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Abstract 

^ ■ 

Q_i' This paper is devoted to the special case of the restricted circular three-body 

^ . problem, when the two primaries are of equal mass, while the third body of neg- 

%^ I ligible mass performs oscillations along a straight line perpendicular to the plane 

c/3 ' of the primaries (so called periodic vertical motions). The main goal of the paper 

^ is to study the stability of these periodic motions in the linear approximation. 

A special attention is given to the alternation of stability and instability within 
^ I the family of periodic vertical motions, whenever their amplitude is varied in a 

' continuous monotone manner. 



1 Introduction 



The term "Sitnikov problem" appeared originally in the context of studies of oscilla- 
tory solutions in the restricted three body problem. These studies were initiated by 
Sitnikov [19]; they stimulated the application of symbolic dynamics in celestial me- 
^ ! chanics [1]. We recall that Sitnikov considered the case when two primaries have equal 

I masses and rotate around their barycenter O , while the infinitesimal third body moves 

along a straight line normal to the plane defined by the motion of the primaries and 
passing through O (usually the motions of the third body perpendicularly to the plane 
of the primaries are called "vertical"; below we will follow this tradition). 

Sitnikov concentrated his attention on phenomena taking place when the primaries 
move in elliptic orbits. More bibliography on "elliptic" Sitnikov problem can be found, 
for example, in [9l [TOl [I2] . 

If the primaries move in circular orbits, then the vertical motions are integrable. 
The corresponding quadratures were presented at the beginning of the XX century 
by Pavanini [T6] and MacMillan [13] - much before the start of Sitnikov's studies. 
Relatively simple formulae for the vertical motions, written in terms of Jacobi elliptic 
functions, can be found in [1]. 

Since the integrability of third body motion is something extraordinary within the 
restricted three body problem, many specialists investigated the properties of vertical 
motions in the case of primaries moving on circular orbit. Very often the term "circular 
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Sitnikov problem" is applied to describe this field of research. Taking into account its 
popularity, we will use it too. Nevertheless, some authors prefer terms like "Pavanini 
problem" or "MacMillan problem" , which are probably more correct from the historical 
point of view. 

Depending on the initial values, three types of vertical motions are possible in the 
circular Sitnikov problem: the hyperbolic escape (i.e., the escape of the third body 
with non-zero velocity at infinity), the parabolic escape (i.e., the escape of the third 
body with zero velocity as the limit at infinity) and, finally, the periodic motion, in 
which third body goes away up to a distance a from the plane defined by primaries 
and then returns to it. 

The first stability analysis of the periodic vertical motions in the circular Sitnikov 
problem was undertaken by Perdios and Markellos [18], but they drew the wrong con- 
clusion that vertical motions are always unstable (Perdios and Markellos only analyzed 
the vertical motions with the initial conditions such that a < 4 ; as it was estab- 
lished lately it is not enough to put any hypothesis about the stability properties of 
the motions with larger values of a). The mistake was pointed out in [4], where the 
alternation of stability and instability of vertical motions were found numerically in 
the case of continuous monotone variation of their amplitude a . Lately the existence 
of such an alternation was confirmed by the results of computations presented in [17] 
and [2D]. Taking into account their numerical results, the authors of [20] proposed the 
hypothesis that the lengths of stability and instability intervals have finite limits as a 
increases. This hypothesis was formulated on the basis of computations in which a did 
not exceed the value 13 . Our numerical investigations demonstrate that the rapidly 
decreasing difference of the stability intervals at a ~ 13 is a manifestation of a local 
maximum of their lengths; if a is increased further, then the lengths of the stability 
and instability intervals tend to zero. 

There is one more important property of vertical motions, which can be observed 
only for a ^ 1 : the intervals of "complex saddle" instability, when all eigenvalues of 
the monodromy matrix are complex and do not lie on the unit circle. According to 
our computations first such an interval begins at a 546.02624 , its length is ~ 10~^ . 
It means the erroneous of the statement in [4] (p. 113), that the stability indexes of 
the vertical motions in circular Sitnikov problem are always real (this statement was 
based on the results of numerical studies in which the amplitude of the motion a was 
smaller 17; as one can see it was not enough for such a general conclusion). 

To conclude our short review on previous investigations of vertical motions' sta- 
bility in circular Sitnikov problem we would like to mention the generalization of this 
problem for systems of four and more bodies [31 EI] . Numerical results presented in 
[3t |2T] demonstrate that in the generalized problem the absence of stability/instability 
alternation in the family of vertical motions persists. 

The aim of our paper is to study the stability property of the periodic vertical 
motions at large values of the "oscillation amplitude" a , both numerically and ana- 
lytically. A special attention will be given to the phenomenon of infinite alternation of 
stability and instability in this family. 

In fact, the infinite alternation of stability and instability in the one-parameter 
family of periodic solutions is rather typical for Hamiltonian systems, although the 
general investigation was carried out only for 2D0F systems [HI E] . Different examples 
can be found in [3 [TTl [T^ . 

Nevertheless, an important difference exists between the circular Sitnikov problem 
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and other systems in which the alternation of stabihty and instabihty was estabhshed 
earher. In the circular Sitnikov problem the discussed family of periodic solutions 
possesses as a limit unbounded aperiodic motions - parabolic escapes, while in pre- 
viously considered systems the corresponding families and their aperiodic limits were 
bounded [3 HB] . Due to this difference, the alternation of stability and instability in the 
circular Sitnikov problem can not be studied in the same way as it was done in [71 [TTl [T3] 
(one could try to compactify the phase space by means of certain changes of variables, 
but we were unable to find any reduction to what was investigated already). 

This paper is organized as follows. In Sect. 2 some general properties of the vertical 
motions are discussed. In Sect. 3 we present the linearized motion equations used in 
our studies of the vertical motions' stability. The results of the numerical investigation 
of the stability are reported in Sect. 4. In Sect. 5 we prepare for the analytical 
investigation: the approximate expression for the monodromy matrix is derived here. 
Using this expression, some important stability properties of vertical periodic solutions 
with large amplitudes a are established in Sect. 6 (in particular, the asymptotic 
formulae for the intervals of stability and instability are obtained). In Sect. 7 we 
discuss briefly the vertical motions in the generalized circular Sitnikov problem with 
four and more bodies. Some concluding remarks can be found in Sect. 8. 



2 Preliminary. Some general properties of the ver- 
tical motions in the circular Sitnikov problem 

We consider the restricted, circular, three-body problem with primaries having equal 
masses, say mi = m2 = m . Let OxiX2X^ be a synodic (rotating) reference frame 
with the origin at the barycenter O ; the masses mi and m2 are arranged on the 
axis Oxi , while the axis Ox^ is directed along the rotation axis of the system. The 
coordinates of the infinitesimal third body in the synodic reference frame will be used 
as generalized variables: 

qi = xi, q2 = X2, qs = X3. 

Below we assume that all variables are dimensionless. 

The equations of motion of the third body can be written in Hamiltonian form with 
Hamiltonian function |1] 

'H = ^{pI+pI+ pI) + Piq2 - P2qi 

Here ri and r2 denote the distance between the third body and the corresponding 
primary, while Pi,P2,P3 are the momenta conjugated to qi,q2,q3 ■ 
The phase space V = {{p,q)} possesses a manifold 

V = {(p,g),Pi = P2 = qi = q2 = 0} 

which is invariant with respect to the phase flow. The phase trajectories lying on V 
correspond to vertical motions with the third body staing always on the axis Ox^ . 
Consequently, the vertical motions are governed by a reduced IDOF system with 
Hamiltonian 

n = 4-^. (1) 
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Figure 1: Phase flow on the manifold V . Thick lines denote the separatrices ("H = ). 

The phase portrait of the system with the Hamiltonian (1) is shown in Fig. 1. It is 
remarkable that the separatrices (the borders between trajectories representing periodic 
motions and hyperbolic escapes) intersect at infinity. 

The periodic solutions associated to the system with Hamiltonian l-L form a one- 
parameter family 

P3{t,a),qs{t,a), (2) 

where as parameter a one can choose the "amplitude" of the periodic motion (i.e., 
a = maxfgfli | gs] ) or the absolute value of at the passage trough the barycenter O 
or the value h of the Hamiltonian "H in this periodic motion. The first variant is the 
most convenient for us, therefore a in (2) will denote the "amplitude" of the periodic 
motion. For definiteness we assume that 



P3(0,a) = 0, ^3(0, a) = a. 

There exist explicit expressions for the solutions (2) in terms of Jacobi elliptic 
functions Since they are not used in the forthcoming analysis, we do not rewrite 
them here, except for the formula about the period of vertical motion: 



T 



V2 



1 - 2A;2 



E{k) + 



TT 



2J2(1-2P) 



1 — Aq l^arcsin 



'1-2P 
1-P ^ 



(3) 



Here E{k) is the complete elliptic integral of the second kind, Ao(v5, k) is the Heuman 
Lambda Function, while the value of the modulus k is given by the formula 



k 



2 + h, 



where 



For motions with large amplitudes ( a ^ 1 ) the following approximate formula can 
be used in place of (3): 

T ^ y27ra^/l (4) 
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As it was mentioned before, the separatrices = {{pf (t) , qf (t)) , t G R^} , repre- 
senting the parabohc escapes, can be interpreted as a formal hmit for periodic motions 
at a — 7- oo . The parabohc escapes obey the approximate law 



qt{t)^±{^\ e'\ (5) 




Formulae (4) and (5) are easily obtained if one suitably relates the properties of vertical 
motions with the properties of rectilinear motions of a particle in a Newtonian field. 



3 The stability problem for periodic vertical mo- 
tions 



Our efforts are concentrated on the analysis of the vertical motions' stability with 
respect to "horizontal" perturbations, due to which the third body leaves the axis 
Ox^. Under the linear approximation, the behavior of the variables pi,P2,qi,q2 in 
the perturbed motion is described by the linear Hamiltonian system of equations with 
periodic coefficients: 

JH(t)z. 



dt 



(6) 



Here 



(Pi,P2,qi,q2] 




H(t) 



/ 1 
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The symbol E^. is used to denote the identity matrix of the k -th order. The function 

1 /2 

D{t,a) = (^ql{t,a) + depends periodically on time with a period = , 

where T{a) denotes the period of the particular vertical motion whose stability is 
investigated. 

As it is known, the restricted circular three-body problem admits several types of 
symmetry (for example, they are used for the numerical construction of 3D periodic 
solutions [IT]). The consequence of these symmetries is the following property of the 
variational equations (6): if z{t) is a solution of (6), then these equations admit the 
solution 

z(t) = Qz(-t), (7) 

where Q it the 4 x 4-diagonal matrix, Q = diag(l, —1, —1, 1) . 

According to Floquet theory, in order to draw a conclusion about the stability or 
instability of the solutions of (6), one should analyze the spectral properties of the 
monodromy matrix M = W(T*,0) , where W(t, t') denotes the normal fundamental 
matrix corresponding to the system (6) (i.e., the matrix solution of (6) with the initial 
condition W{t',t') = E4). 

The normal fundamental matrix corresponding to the linear Hamiltonian system 
(6) is a symplectic one, i.e. 



W^(t,t')JW(t,t') = J. 
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It is also worthwhile to mention some other properties of this matrix: 

W(t, t") = W{t, t')W{t', t"), W{t + T„t' + n) = W(t, t'), (8) 

W(0, -t) = QW-\t, 0)Q = -QJW^(t, 0)JQ. 

The first two equalities in (8) are elementary, while the last one is a consequence of 
the symmetry property (7). 

Using the relation (8) one easily obtains 

M = QW-i o) QW o) = -QJW^ o) JQW o) . 

The characteristic equation of the system (6) 

det (M - PE4) = (9) 
is reciprocal and it can be written as 

- cip^ + C2p^ - cip + 1 = 

where 

3 4 

Ci = trM , C2 = ^ ^ {nijjmkk - mjumuj) . 
j=i k=j+i 

The quantities m^j in the last formula are the elements of the monodromy matrix M . 
It is also possible to rewrite the characteristic equation (9) as the product 

(p'-26ip + l)(p2_262P + l) = 0. (10) 

The coefficients 61, 62 h^ (10) are the roots (real or complex) of the quadratic equation: 

Ax^ - 2cix + (C2 - 2) = 0. 

Often enough the quantities 61,62 are called the stability indices [3]. The periodic 
vertical motion is stable whenever {61,62} C / = (—1, 1) C (i.e., when 61,62 are 
real and their absolute values are smaller than 1). In the case 

{61,62} C/= [-1,1], {6i,62}(// 

an additional investigation is needed to draw a conclusion about stability or instability. 
In all other cases the instability takes place. 

4 Numerical results 

We recall that in f4j the alternation of the stability and instability in the family of 
periodic vertical motions (2) was discovered. Later on, more accurate results were 
published in pO]: the length of the first 35 intervals of stability and of the first 34 
intervals of instability was calculated. In |2U] also an attempt was undertaken to 
establish certain regularity in the variation of these quantities: the existence of non- 
zero limits for the intervals' lengths was proposed. 

In Fig. 2 and Fig.3 we present the results of some calculations, when the first 700 
intervals of stability and instability are considered. The graph in Fig. 2 shows that for 
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Figure 2: Length of the stabihty interval as a function of its number 



the first 30 intervals of stability the length of the intervals increases and only afterwards 
the decrease of the length takes place. The hypothesis formulated in [20] was based on 
the wrong interpretation of the small variation of the intervals length in vicinity of the 
maximum. In Fig. 3 the length of the instability intervals decreases monotonically and 
it does not follow the empirical law derived in [20] (according to this law, the length 
of the instability intervals has the limit Aj„st ~ 0.254 ; evidently, it is not so). 

Our results allow us to propose the following approximate formulae to characterize 
the behavior of the stability and instability intervals' length in Fig. 2 and Fig. 3: 

Ast-0.25N-'/\ Ai„,t ^ 0.584iV-i/3. (11) 

More precisely, these formulae are valid for the periodic vertical motions with amplitude 
a smaller the critical value = 546.02624... The reason of such a restriction and the 
situation for a > a* will be revealed a little bit later. 

It is also useful to discuss here in what way the length of the stability intervals 
Ast and the length of the instability intervals Ainst depend on the amplitude of the 
vertical oscillations. Under the same restriction a < a* we obtain from our numerical 
investigations 

A,t ~ 0.3a-'/\ Ainst ~ QMa-^'\ 

Remark. If one needs a rigorous definition about the meaning of the quantity a in 
the last formulae, one could interpret it as the boundary value between two successive 
intervals of stability and instability. 

In Fig. 4 the behavior of the coefficients 61,62 appearing in the characteristic 
equations (10) is shown. Fig. 4a, 4b and 4c allow us to compare the properties of these 
coefficients, when the parameter a varies in different intervals. All graphs demonstrate 
the approximate periodicity, their period with respect to the parameter a corresponds 
to an increase of period of vertical oscillations T of about Svr . It is important to point 
out the small gaps in the Fig. 4c: for the corresponding value of the parameter a (i.e., 
when a belongs to the intervals where the graphs are not defined) the stability indices 
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Figure 3: Length of the instabihty interval as a function of its number N (the first 
interval is not presented: if it was shown in the same scale with all subsequent intervals, 
it would have been difficult to understand the behaviour of the graph for large N ) 

6i, 62 have complex values and the so-called "complex saddle" instability of the vertical 
motion takes place. The enlarged fragments of the graphs in the vicinity of the gaps 
are given in Fig. 5. 

As it follows from our calculations, the first interval of " complex saddle" instability 
begins at a = a* . Since such a value of vertical motion amplitude is large enough, it 
provides us with an explanation why this kind of instability of vertical motions in the 
circular Sitnikov problem was not recognized in previous studies where relatively small 
values of a were considered. 

Increasing further the parameter a (i.e., for a > a*), we observe a stabil- 
ity/instability alternation of more complicated type: "wide" interval of instability - 
"narrow" interval of stability - "narrow" interval of "complex saddle" instability - 
"wide" interval of stabihty - "wide" interval of instability - ... An analog of the for- 
mulae (11) can be constructed in the case a > a* , but we prefer to present in Sect. 6 
several asymptotics written in a more convenient way. 

Finally it worth while to mention that the Runge-Kutta-Fehlberg method of 7-8 
order with variable step was used to integrate numerically the variation equations (6). 
The accuracy of the integration procedure (the local tolerance) was taken 10~^° . Since 
the period of vertical oscillations increases proportionally a^^^ the variation equations 
should be integrated over relatively large time intervals: if we take for example a = 500 
then the value of half-period ^ 2.4837 ■ 10"^ . To check the influence of the round- 
off errors some computations were done both with double and quadruple precision 
arithmetic. 
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Figure 4: The behaviour of the coefScients bi and 62 appearing in the characteristic 
equation (10) 
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Figure 5: The enlarged fragments of the coefficients graphs (see Fig. 4) in the vicinity 
of the gaps 
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5 Approximate expression for monodromy matrix 



In this section an approximate expression for the monodromy matrix M is derived. It 
will be used to discuss the phenomena described in Sec. 4 (the alternation of stability 
and instability, the decrease of stability and instability intervals by increasing the 
parameter a , etc) . 

We assume that the amplitude a of the periodic solution (2) is so large, that we 
can define an auxiliary quantity d such that 



1 (i < a. 



(12) 



To start with we write down the monodromy matrix M = W(T*, 0) as the product 
of three fundamental matrices: 



M = w(r.,t,-)w(t,-,t+)w(t+,o), 



(13) 



where tt G 



0,f 



and r„ 



— are the instants at which the third body is at 



distance d from the barycenter O in the periodic vertical motion (2) (at t — the 
third body moves away from barycenter, at t — it approaches the barycenter). 

Approximate expression for the matrix W(tJ, 0) . If the condition (12) is satisfied 
the phase point (psit) , q^it)) moves on the manifold V in close vicinity of the scparatrix 
S*"^ at t G [0,t^] . Within such time interval the difference between q^{t^a) and qt{t) 
is small enough. Neglecting this difference, we replace q^it^a) in (6) by qtif) ; as a 
consequence, the normal matrix solution W+(t, 0) of the obtained system provides us 
the suitable approximation for W(t, 0) at t G [0, f^] . 

The behavior of W+(t, 0) at t — > +oo is described by the remarkable asymptotic 
formula: 

W+(t,0)ftiR(t)A(4(t))U. (14) 



Here 



R(t) 



( cos t 
— sin t 





sin t 
cos t 










cos t sin t 

■ sin t cos t J 



I 



U 



A(g3) 



/ 0.3248 . . 
0.1302.. 
1.1175.. 
V 0.2113.. 



1 

93 



Tf # 







V 



1 

93 




^ 


0.1020, 
0.1189, 
0.1718, 
0.6404 , 





-93 


-0.4664. 
0.5296, 
1.4408, 
0.9414, 



-93 / 

0.2228, 
-2.0211. 

0.4791 . 
-2.5646 . 



The derivation of the formula (14) is based on some simple ideas. Let us take o? 1 
and write down W+(i, 0) as the product 



W+(t,0)=W+(t,%)W+(t^,0), 



(15) 
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where t-^ is the moment of time when the third body is at distance d from the barycen- 
ter O in the motion corresponding to the parabohc escape — qtit) ■ As next step, 
we modify the equations (6) to find the approximate expression for W+(i, %) at t > t-^- 
Since at t > the third body is far enough from the primaries mi and m2 , it looks 
natural to replace D by q^{t) in the right parts of the first two equations in system 



(6) and to neglect the small term 



dz 

'di ' 



The system (6) takes the form 



JH(g3+(t))z 



(16), 



with 



H(g3) 



/ 1 








1 -1 
-1 







1 





1 






Now it is worthwhile to make the following remark. Let us consider the rectilinear 
parabolic escape of the material point in the field of an attracting center. Under a 
proper choice of units, the distance between the attracting center and the point varies 
as 



(17) 



If the asymptotics (5) is used for q^{t) in the equations (16), then these equations 
coincide with the motion equations of the above mentioned material point, linearized 
in the vicinity of the solution (17) and written in the reference frame uniformly rotating 
around the line of the escape. 

Taking this into account, we implement in (16) the change of variables 



where 



Z = (Pl,P2,gi,g2) Z = (pi,P2,?l,?2) 



z = R(t^-i)z. 



This change of variables can be interpreted as the transfer from the synodic reference 
frame 0x^X2X3 to the sidereal (fixed) reference frame OxiX2X^ (OxsHOxs). As a 
result the linearized equations of motion split into two independent subsystems 



dt 



— ' 



dqi _ 



i = 1,2. 



(18) 



It is not difficult to find partial solutions to the system (18) 



Pi 




Qi = ?3, P3- 



0, Qs-i = 0, 



1,2 



and 



1 

Pi = — , 

?3 



Qi = QsQs 



'2^3, Ps-i = 0, qs-i = 0, 



1,2. 



Here and below the dots are used for derivatives with respect to time. 
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Four independent partial solutions allow us to write down the normal fundamental 
matrix in terms of the variables z : 

Coming back to the initial variables, we get 

W(t,t^) = R(t-t^)W+(t,t^). (19) 

Substituting (19) into (15) we obtain the expression for the normal fundamental 
matrix W+(t, 0) as the product of three matrices with only one of them depending on 
time: 

W+(t,0)^R(t)A(g3+(t))U(d). (20) 

Here 

U(d) = A-i(S)R(-t^)W(t^,0). 

The formula (20) can be used to compute the elements of the matrix W_|_(t,0) at 
t 1 . Asymptotically their values should not depend on the choice of d. It means 
that the following limit exists: 

U = _lim \j(d). 

Substituting U instead of V{d) into (20) we arrive to the formula (14). 

The fundamental matrix W_|_(t, 0) was introduced in such a way that it provides 
the vertical motions satisfying (12) with a "universal" (i.e., independent on a) ap- 
proximation W(t, 0) ~ W+(t, 0) at t G [0,t^]. Using the relation (14), we finally 
obtain 

W(t+,0)^R(t+)A(rf)U. (21) 

Approximate expression for the matrix W(t^,tj") . Since at t G [t^j^d] ^^e third 
body is far enough from the primaries, we neglect again the difference between their 
gravity field and the gravity field of the attracting center placed at the barycenter 
O . To obtain the expression for W(t^,tj|) within such an approximation we need to 
integrate the system 

^ = JH(g3(t,a))z, (22) 

where q^{t^a) describes the motion in the Newtonian field along the segment [0,a] 
on the axis Ox^ . It is supposed that the maximum distance a from the body to the 
attracting center is achieved at t = ^ . In this case qs(t,a) ^ q^{t^a) at t G [t^,^^^] ■ 
Of course the motion along a segment corresponds to the singular impact orbit [22] . 
but it is used here to approximate the regular vertical motion on the time interval were 
the singularities are absent. 
The change of variables 



z = {Pi,P2,qi,q2f ^^ z = {pi,P2,qi,Q2 

where 

z = R/^ 

2 



z = R(^-t)z, (23) 
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allows us to rewrite the equations (22) in the more simple form: 

i = l,2. (24) 
It is easy to check that the system (24) admits the following partial solutions: 



dt 



Qi^ dqi _ 
qi dt 



Pi = 1i = 53, P3- 



0, q3-i = 0, 1 = 1,2 



(25) 



and 



1 2 ^ _ 

Pi^- , Qi^ qsqs: Pa- 

qs a 



0, g3_, = 0, i = l,2. 



(26) 



To compute q^ in (25) and (26) the energy integral can be used. In the case of the 
motion along the segment [0, a] in the Newtonian field, this integral takes the form 



^2 

2 



1 

?3 



1 

a 



and consequently 



Taking into account (25), (26) we write down the fundamental matrix for the system 
(24) as 



W U. 



T 















2 - 



Q3 



-aqsq^ 
\ -a?3?3 



^ 

?3 



^ 



23 

a 



and then (taking into account the relation (23)) we write the matrix for the system 
(22) 



T 



Using the expression (27) we find 



W (t, — ) =R(t- — ]W it,— 



t,,':^)=R(t,-'^)N{d,a), 



(27) 



W 



where we denote by 



2 ) 



N(d,a)=W(t^,^ 



-If^) 







^2 - 1) 



4./2 



1 _ 1 



a 





d 
a 
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The final step is based on the last formula in (8), namely 



Here 



K((i) 



/ -3 2^d-3/2 \ 

-3 2v^(i-3/2 

2^2^3/2 0-3 

V 2V^(i3/2 -3 y 

For completeness we should add the following formula: 



K(rf) ^ N(d,a)QN-^(d,a)Q. 

Approximate expression for the matrix 'W{T^,t'2) ■ Using again the relations (8) 
we get 

W(T.,t,-) = W(0,-t+) = QW-^(t+,0)Q. (29) 

Then the substitution of the previously obtained expression for W(t^, 0) (the formula 
(21)) into the right part of (29) provides us with the desired approximate formula for 
W(T„t,-). 

Finalizing the construction of the approximate formula for the monodromy ma- 
trix M. The substitution of the approximate expressions for W(T*,i^) , W(i^,iJ) , 
W(i 0) into (13) yields 

M(a) ^ {Q (A(ci)U)-^ QR(T, - t,")} {R(t, - t+)K(ci)} {R(t+)A(rf)u} . (30) 

Simplifying (30) the extraordinary simple result can be obtained 

M(a) ^ QU-^QR(T,(a))U. (31) 

This formula allow us to investigate analytically the stability properties of the periodic 
vertical motions in the case a':^ 1 . 



6 New insight into the stabihty properties of the 
vertical motions 

As it follows from (31) the coefficients of the monodromy matrix M and, respec- 
tively, the coefficients of the characteristic equation (9) are 27r -periodic functions of 
the semiperiod of the vertical motion T^, . To illustrate this we present in Fig.6 the 
graphs of the coefficients 61,62 (only the real values) as varies in the interval 
[27rn, 27r(n + 1)] , where n is a large enough integer number. Using the formula (3), 
which defines the dependence of on the amplitude a , it is not difficult to prove 
that in terms of a the lengths of the stability and instability intervals decrease pro- 
portionally to a~^/^ as a — > -|-oo . 

Taking into account the approximate expression for the monodromy matrix M , 
we describe in more details the repeating pattern of stable and unstable intervals 
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Figure 6: Graphs of the coefficients 61,62 computed on the base of the approximate 
formula for the monodromy matrix M . Only real values are shown. 

mentioned at the end of Sec. 4. This pattern consists of four intervals appearing 
in the following order as a increases: 

"Wide" interval of instability. Both coefficients 61, 62 are real, but one of them has 
absolute value greater than 1 ("saddle-center" instability). The asymptotic length of 
the interval in terms of the amplitude of the motion is about 0.643544 • a~^/^ , while 
the variation of the semiperiod T* equals to about 2.144392 . 

"Narrow" interval of stability. The coefficients 61,62 are real and belong to the 
interval (—1, 1) . The approximate length is 0.068655 • a~^/^ ; the variation of the 
semiperiod equals to about 0.228768 . 

"Narrow" interval of instability. The coefficients 61,62 are complex ("complex 
saddle" instability). The approximate length is 0.048166 • a~^/^ ; the variation of the 
semiperiod equals to about 0.160497. 

"Wide" interval of stability. The coefficients 61,62 are real and belong to the 
interval (—1,1) again. The approximate length is 0.182445 • a"^/"^ ] the variation of 
the semiperiod equals to about 0.607936 . 

To conclude, we recall that before the ffist appearance of the interval of "complex 
saddle" instability at a = a* , a more simple pattern with only two intervals was 
observed. The "transient" asymptotics for the length of the stability intervals in the 
case 1 <^ a < can be obtained by adding of the lengths of the "narrow" instability 
interval and both stability intervals in the final pattern. It yields 



which is in good agreement with the corresponding numerical result presented in Sec. 
4. 

7 Stability of the vertical motions in the circular 
Sitnikov problem with four and more bodies 

The investigation of the generalized circular Sitnikov problem with four and more 
bodies revealed that in contrast to the case of the three body problem there is no 
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alternation of stability/instability in the family of vertical motions [3| [2T]. 

For simplicity we limit our consideration to the case of the restricted four body 
problem. It is assumed that three primaries of equal mass rotate around the barycenter 
O in circular orbit with the radius R = l/\/3 [21]. Under the linear approximation 
the stability analysis of the fourth body periodic vertical motion q{t, a) is reduced to 
the study of the spectral properties of the monodromy matrix associated to the system 
of linear differential equations with periodic coefficients 



dz ^ , ^ 
- = JHWz. 



(32) 



where 
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D{t,a) = [ql{t,a) + - 
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It is remarkable that equations (32) possess a circular symmetry: for any real a 
they are invariant with respect to transformations of the form 

z = Rfa)z, 



while the non-linearized equations of motion of the fourth body in the synodic reference 
frame admit only the rotational symmetry of the 3rd order. The possibility for the 
linearized equations of motion to have a larger group of symmetries in comparison 
to the original non-linear system was pointed out by V.I. Arnold [2] (Sec. 23). In 
particular, in the case of the initial system rotational symmetry of -th order ( iV > 3 ) 
the linearized equations always have circular symmetry. This is the reason why the 
stability analysis of the vertical motions, based on the linearized equations, yields 
similar results for the Sitnikov problem with four and more bodies and for the particle 
dynamics in the gravity field of the circular ring [5] (numerically it was shown in [3J). 

It is convenient to rewrite the equations of motion (32) in a sidereal (fixed) reference 
frame by means of the transformation of variables 



{Pi,P2,qi,q2) 



z 



where 

z = R(-t)z. 

After that the equations of motion split into two identical independent subsystems: 



dpi 
dt 



qi 



2D3 



dqi 



1,2. 



(33) 



Let W_(_(t, t') denote the normal fundamental matrix for the system (33) in the 
case when ^it.a) is replaced by q^(t) , which corresponds to the parabolic escape. 
Using the same technique as in Sect. 5 we obtain the asymptotic formula 



W+(t,0)^A(g3+(t))U, 
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where 



1.2690 \ 

-1.2690 
0.7061 
-0.7061 / 

Applying the main ideas of Sect. 5, we estabhsh the following property of the 
monodromy matrix M(a) associated to (33): at a — )■ +oo the matrix M(a) M* , 
where the constant matrix M^, = QU^^QU. The eigenvalues of the matrix M^, are 
the asymptotic limits for multiplicators (of multiplicity 2) of the system (33): 

lim Pi (a) = p*, p* = -0.4446 . . . , pX = -2.2488 . . . 

Finally, it is not difficult to derive the asymptotic formulae the for multiplicators 
of the original system (32): 

pi,2 ~ Piexp(±zTj, p3,4 ^ p;exp(±iT,). 

On the complex plane pi, . . . ,p4 are placed in the small vicinity of the circles with 
radii |pi| < 1 and |p2| > 1 . Consequently in the circular Sitnikov problem with four 
bodies the periodic vertical motions with large amplitudes are always unstable. 

Finally we would like to note the another opportunity to introduce the generalized 
circular Sitnikov problem with N bodies using the appropriate straight line solution 
of the problem of (A^ — 1) bodies |T3]. If in such a solution (A^ — 1) primaries are 
arranged symmetrically with respect to the barycenter then the infinitesimal A^ th body 
can move periodically along an axis around which the rotation of the primaries takes 
place (naturally, the proposed generalization is possible for odd A^ only). Likely this 
family of periodic motions exhibits the alternation of stability and instability. 



U= lim A-\d)W+{td,0) 

d—> + OD 



/ 0.2456 

0.2456 

0.9246 

V 0.9246 



8 Conclusion 

The combination of numerical and analytical approaches provided us with the opportu- 
nity to correct, clarify and extend some previously known results related to the circular 
Sitnikov problem (mainly about the stability of vertical motions). For the first time 
under the scope of this problem the possibility of the "complex saddle" instability was 
revealed within the family of vertical motions. 

For our theoretical constructions it was essential that the phase trajectories corre- 
sponding to the solution under consideration have lengthy parts in the vicinity of the 
peculiar separatrices of the problem - the parabolic escapes to infinity. Often enough 
it is possible to introduce a suitable auxiliary mapping in the vicinity of the separatrix 
in order to study the local properties of the phase fiow. It would be very interesting to 
develop similar for the circular Sitnikov problem. 
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